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Charging of a clean two-dimensional island is studied in the regime of small concentration of 
electrons when they form the Wigner crystal. The number of electrons in the island is assumed to 
be not too big (N < 100). It is shown that the total energy of the island as a function of N has a 
quasi-periodic component of a universal shape, that is independent of the form of electron-electron 
interactions. These oscillations are caused by the combination of the geometric effects associated 
with packing of the triangular lattice into the circular island. These effects are: the shell effect, 
associated with starting a new crystalline row, and the so-called confinement polaronic effect. In 
the presence of close metallic gates, which eliminate the long-range part of the electron-electron 
interactions, the oscillations of the energy bring about simultaneous entering of the dot by a few 
electrons. 

PACS numbers: 73.20.Dx, 73.40. Gk, 73.40. Sx 
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I. INTRODUCTION 

In recent experimcntsBi the charging of a quantum dot 
is studied by the single electron capacitance spectroscopy 
method. The quantum dot is located between two capac- 
itor plates: metallic gate and heavily doped GaAs layer. 
Tunneling between the dot and the heavily doped side 
is possible during the experimental times while barrier 
to the metal is completely insulating. DC potential V g 
and a weak AC potential are applied to the capacitor. 
With the increase of V g the differential capacitance ex- 
periences periodic peaks when addition of a new electron 
to the dot becomes possible. The spacing between two 
nearest peaks AV g can be related to the ground state 
energy E(N) of the dot with N electrons: 



aeAV g = E(N + 1) - 2E(N) + E(N - 1) 



A(N) = e 2 /C, 



(1.1) 



AT • 



Here a is a geometrical coefficient, A(iV) is the charging 
energy, Cm is the capacitance of the dot with N electrons. 
It was observed in Refs. |^J| that at a low concentration of 
electrons or in a strong magnetic field the nearest peaks 
can merge, indicating that at some values of V g two elec- 
trons enter the dot simultaneously. In other words some 
charging energies apparently become zero or negative. In 
a fixed magnetic field this puzzling event is repeated pe- 
riodically in N. Disappearance of the charging energy 
looks like a result of unknown attraction between elec- 
trons and represents a real challenge for theory. 

Explanation of the pairing of the differential capac- 
itance peaks based on the lattice polaronic mechanism 
has been previously suggested in Ref. ||. In Ref. || it was 
demonstrated how electron-electron repulsion screened 
by a close metallic gate can lead to electron pairing for a 
specially arranged compact clusters of localized states in 
a disordered dot. This effect is a result of redistribution 



of the other electrons after arrival of new ones. It was 
interpreted in Ref. [I| as electronic bipolaron. 

In this paper we study the addition spectrum of a dot 
where density of electrons is small and disorder is weak 
enough for electrons to form the Wigner crystal. We 
call such a dot the Wigner crystal island. In the experi- 
mental conditions of Ref. [l] one can think about Wigner 
crystal island literally only in the highest magnetic field. 
One can also imagine similar experiments with a Wigner 
crystal island on the surface of liquid helium. 

We show that when a metallic gate is very close to the 
plane of a two-dimensional (2D) island, two or more ad- 
dition levels can be anomalously close and even merge 
similar to the observations of Ref. [I]. This phenomenon 
repeats periodically, and is associated with the formation 
of highly symmetric electronic configurations, consisting 
of an integer number of closed radial shells. The new elec- 
trons added to such a configuration reside on its surface. 
They exert some force onto the symmetric configuration 
making it move in the external confinement field. This 
displacement facilitates the entrance of the new electrons 
mediating an effective attraction between them. To em- 
phasize that this effect is produced by the confinement 
potential and to stress its multi-electron nature we call 
it confinement multi-polaronic effect. 

The paper is organized as follows. In Section || we 
present our model and the resu lts of its numerical solu- 
tion. In Sections III and IV we give the explanations 
of the numerical results for the screened and for the un- 
screened Coulomb interactions respectively. Section |v| 
is dedicated to our conclusions. 



II. NUMERICAL SIMULATIONS 

In this Section we consider the system of N interact- 
ing electrons in the parabolic confinement. The energy 
of the system is given by the following expression: 



1 



E = Y J U{r l 



+ A 



(2.1) 



Here A plays the role of strength of the confinement, 
and U (r) is the interaction potential. The forms of the 
interaction potential considered are: the pure Coulomb 
interaction 



U (r) = e 2 /K,r, 



(2.2) 



with e and k being the electron charge and the dielectric 
constant correspondingly, and the exponential interac- 
tion 



U (r) = U Q exp {-r/d) . 



(2.3) 



The latter interaction potential corresponds to the case 
when the island is situated between two metallic gates, 
with nd being the distance between them and Uq ~ 
e 2 jnd. 

This energy was minimized over the coordinates of 
electrons using the numerical technique similar to that 
outlined in Ref. |[ It is the standard Metropolis simu- 
lated annealing algorithm. The temperature of the sim- 
ulation was decreased exponentially (quenched) in 10 5 iV 
simulation steps. The average displacement of electrons 
in a step was chosen automatically to assure the accep- 
tance probability of a new configuration of 0.3. Each 
time the simulation for the new number of electrons was 
started from a random set of coordinates. We used a few 
restarts for a given number of electrons. Although the 
configurations obtained in different restarts sometimes 
were different, the maximum difference in energy is typ- 
ically 10 times smaller than the systematic variations of 
E(N). Thus although we cannot guarantee that we have 
found the ground state of the electron system, we are con- 
fident that these variation of E (N) are reproduced reli- 
ably. Our results for energies agree with those of Ref. |5| 
and for some N are even lower. Before presenting these 
results we would like to discuss the method of the data 
processing. 

Consider first the problem in the case of n o sc reen- 
i ng, that is described in our model by Eqs. (2.1) and 
(|2.2[) . The instructiveness of this problem is in the fact 
that the corresponding electrostatic problem (ignoring 
the discreteness of the charge of electron) can be solved 
exactlyB. The solution for the density of electrons can be 
shown to be hemisphere: 



n(r) = n \ 1 



R 2 
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(2.4) 



The total potential of the electrons is equal to tp = 
2AR 2 . This solution allows us to calculate the scaling 
of all other contributions to the energy of the ground 



state with the number of particles (see Appendix |A|) . In 
particular: 



E(N)= (e 2 / K ) 2/3 A 1 / 3 (r ]1 N^ 3 + 

+ m N 7 ' & + r) 3 N 2 / 3 + m N 7 ^ 15 + ...), 



(2.5) 



where 77.; are some constants. The first term in this 
asymptotic series is the electrostatic energy that can be 
found from the solution cited above. The next three term 
in the order of appearance are: the correlation energy, 
the overscreening energy associated with the screening 
of the external potential by the Wigner crystal, and the 
surface energy. The coefficients r\i can be found from 
the best fit to the numerical data: r\\ = 6/5 (37r/8) 2//3 , 
r/2 = -1.0973(2) , 773 = -0.182(3) , 774 = -0.06(1). This 
smooth contribution to energy can now be subtracted 
from the numerical data to obtain the fluctuating part: 



SE (N) = E(N)-E{N). 



The fluctuating part is displayed in Fig. [!]. 



(2.6) 
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FIG. 1. Fluctuating part of energy of the Wigner crys- 
tal island. 

The curve in Fig. |l| evidently has a quasi-periodic 
structure. It consists of the series of interchanging deep 
and shallow minima, separated by the peeks with more or 
less smooth slopes. We call this structure the "Kremlin 
wall" for its resemblance to the embattlements on top of 
the walls of the Moscow Kremlin. Its period with a very 
good precision scales as N 1 / 2 and numerically is equal 
to the number of electrons in the outer crystalline row 
of the island. The amplitude of the oscillations does not 
change with N appreciably. The minima of SE are asso- 
ciated with the highly symmetric electron configurations 
that can be formed at some distinct numbers of electron 
N m that we call magic numbers. One of the examples of 
such configurations N — 85 is displayed in Fig. 0. 
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As one can see, at some numbers of electrons in the is- 
land the total chemical potential is a decreasing function 
oiN. 
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FIG. 2. The magic number configuration N = 85. 

One can use the numerically obtained ground state en- 
ergies to calc ulate the differential charging energy A(iV) 
by Eq. (1.1). The results are shown in Fig. p[ Note 
that although the average charging energy ~ e / kR de- 
creases oc TV" 1 / 3 the relative fluctuations remain w 0.15 
for number of electrons in the well N < 100. 
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FIG. 3. The charging energy of the Wigner crystal 
island as a function of the number of electrons in the ab- 
sence of screening. 

It was surprising for us to find that quasi-periodic os- 
cillations similar to Fig. |l| exist for all the types of in- 
teractions that we have tried. The results for a partic- 
ular ex amp le of the exponential interactions described 
by Eq. 2.3 are presented in Figure || From the to- 
tal energy the smooth component of the form E(N) = 
■qiN 2 + i] 2 N + 7/ 3 iV 2 / 3 + 774/V 1 / 2 + 7757V 1 / 3 is subtracted, 
where the coefficient rji are chosen to minimize the fluc- 
tuations. 

It is clear that the fluctuations retain their quasi- 
periodic form similar to the pure Coulomb case. How- 
ever, contrary to the Coulomb case they grow with N as 
iV 7 , where 7 = 0.8 ±0.1. The relative fluctuations of the 
energy also become bigger because the smooth compo- 
nent of the ground state energy is significantly reduced 
due to screening. As a result at some numbers of elec- 
trons in the island the fluctuations of the charging energy 
exceed the average. The result is demonstrated in Fig- 
ure o showing the total chemical potential of the system 
/j (N) — E (N + 1) — E (N) with exponential interaction. 
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FIG. 4. The fluctuating part of the ground state en- 
ergy of the model with the exponential interaction (2.3) 
for A = W~ 8 ,U = l,d= 1. 

Negative dfi/dN means that two or more electrons "at- 
tract" each other and, with increasing /_t or V g , enter the 
island simultaneously. In other words this means that a 
few Coulomb blockade peaks merge together forming a 
bunch in the charging spectrum. The detailed analysis 
of this instability is given in the next Section. 
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FIG. 5. The total chemical potential n(N) for the 
same parame- 
ters as in Figure | (A = 10" 8 ,i7o = l,d = 1). In the 
vicinity of N = 5, 7, 15, 20, 30, 40, etc. it is a decreasing 
function of number of electrons. 



III. A SIMPLE MODEL FOR THE 
SHORT-RANGE INTERACTIONS 

In the previous Section we saw that bunching is as- 
sociated with screening of the interactions in the island. 
Once these interactions become short-range the fluctu- 
ations of the charging energy exceed the average value. 
To understand the nature of these fluctuations it is in- 
structive therefore to consider the simplest form of the 
short-range interactions: the hard wall interactions. 

Assume that electrons are hard disks of diameter d put 
into the parabolic confinement. The energy of the system 



3 



now is given by 



E 



(3.1) 



where r$ are the coordinates of the centers of the discs. 
This expression formally coincides with the moment of 
inertia of the system of particles of mass A relative to 
the center of the parabolic confinement. Hence, to min- 
imize the energy one has to minimize the correspond- 
ing moment of inertia. As it is well known the moment 
of inertia of a system is minimum relative to the center 
of mass of the system. We conclude that the center of 
parabola is situated in the center of mass of the minimum 
energy configuration. 
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FIG. 6. a) Disk consisting of an integer number of 
closed electron shells and a few electron in the new outer 
shell. One of such inner shells is shown by empty circles. 
New electrons added to the disk (dark circles) are put on 
its surface. In the simplest approach the new electrons 
are spread uniformly over the circumference of the disk, 
b) The fluctuating part of the total energy of the system 
(top) and the total chemical potential (bottom) obtained 
in this approximation. 

Assume that at some number of particles in the well 
they form a highly symmetric configuration consisting of 
an integer number of closed radial shells. These configu- 
rations are observed when the number of particles corre- 
sponds to the deep minima of the fluctuating part of the 
energy of the system (see Fig. [j]). The electrons in these 
magic number configurations form a figure very close to 
the perfect hexagon. If the number of electrons is not too 
big the hexagon is almost indistinguishable from disk. In 
fact the amplitude of the deviations of disk from hexagon 



are less than the radius of particles if the number of par- 
ticles in the well does not exceed 175. Below for the sake 
of instructiveness we replace hexagon with disk. 

When adding a new electron to the symmetric configu- 
ration it has to be placed somewhere on the boundary of 
the island. The moment of inertia of the system relative 
to the center of the disk is then 



E = AN R 2 /2 + A6N(R + dV3/2)< 



(3.2) 



where ANq and R are the mass and the radius of the 
disk, and SN is the number of newly added electrons 
(see Fig. [|). In the first approximation one can think 
that electrons are spread uniformly over the circumfer- 
ence of the disk (Fig. ^a). Then the center of mass of 
the syste m co incides with point O and its energy is given 
by Eq. ([3.2]). It is necessary to notice now that this 
expression is linear in the number of new electrons. At 
the same time the average moment of inertia of the sys- 
tem grows oc iV 2 . Indeed it is equal to NR 2 /2, with R 
being the average radius of the island determined from 
the expression nnR 2 — N, where n — (d 2 \/3/2) is the 



density of centers of particles. Hence the expression (3.2) 
is just a linear approximation to the average energy with 
the difference between them being the fluctuations. The 
fluctuating part of the energy and the total chemical po- 
tential of the system obtained in this approximation are 
shown in Fig. ||b. 

The total chemical potential as a function of the num- 
ber of electrons comprises the series of steps. This de- 
pendence resembles one for the system of non-interacting 
electrons in the conditions of the integer quantum Hall 
effect. The fluctuating part of energy is just a sequence 
of parabolas. The minima of the curve correspond to the 
deep minima of the Kremlin wall structure shown in Fig 

The appearance of the sequence of shallow minima in 
Fig. |l| is associated with displacements of the electrons 
on the surface of the disk. One can notice that these elec- 
trons are actually attracted to each other. Hence instead 
of being uniformly spread over the surface they form a 
compact cluster (see Fig. ^,). This interaction is similar 
to the attraction between two people in a boat. One per- 
son tilts the deck of the boat and another has a tendency 
to slide down the slope. For our model this can be most 
easily understood using the analogy with the moment of 
inertia. Indeed due to the known theorem the moment of 
inertia relative to point O is related to one with respect 
to the center of mass by 
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(3.3) 



Here tcm is the distance from the center of mass to the 
point O. It is clear that by forming a compact group the 
particles on the surface do not change Eo, while at the 
same time they displace the center of mass. But by doing 
this they decrease the total energy, as it follows from the 
previous equation. The results are displayed in Fig. 0. 
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The fluctuating part of energy now looks similar to 
that in Fig. [l]. A simple calculation shows that this 
Kremlin wall structure is described by the following de- 
pendence: 



SE (SN) 
An- 1 



-SN (SN Q - SN) 



4A 

•n-2 



nSN 

~S~Nn 



(3.4) 

where 8 No = 2nRo jd — [v / 2li r /Vo] , is the number of elec- 
trons in the outer shell. By [• • •] we denote the integer 
part of a number. The first term in this equation de- 
scribes the pure shell effect (without attraction), while 
the second is the result of displacement of the center of 
mass (Eq. (O)). 
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FIG. 7. a) Clustering of electrons on one side of the 
disk that shifts the position of the center of mass (point 
CM) from the center of the disk, decreasing the total en- 
ergy (Eq. (3.2)). b) The fluctuating part of the total 
energy (top, bold line) reduced relative to the simplistic 
approach (dashed line) and the total chemical potential 
(bottom, bold line) contrasted to the simplistic approach 
(dashed line). In the regions of negative dfi/dN electrons 
enter the island in bunches. The number of electrons in a 
bunch is determined from the Maxwell's rule, illustrated 
in the inset. 

Another interesting feature of the hard-disk model is 
that the total chemical potential for some fillings of the 
island is decreasing as a function of number of electrons. 
In this case, if the chemical potential in the island is 
fixed by an external gate, the system becomes unstable 
and a few electrons enter the island simultaneously. This 
is a result of the aforementioned attraction between elec- 



trons on the surface of the disk. This attraction can be 
thought of as a multi-polaronic effect based on the spring 
connecting the magic number configuration with the cen- 
ter of the parabolic confinement. Therefore we call this 
effect the confinement multi-polaronic effect. 

Let us examine this instability in more detail. Assume 
that the island is in the proximity of a metallic gate, in 
which the chemical potential is equal to eV g . When AN 
electrons tunnel from the gate to the island the energy 
of the system is changed by the following amount: 

pNo+AN 

AE = n{N)dN -aeVgAN, (3.5) 

JNo 

The first and the second terms in this expression describe 
the change of the total energy of the island and the work 
done by the external source respectively. Introducing 

MAO 

this tunneling is energetically possible: 



fj, (N) — aeVg, we obtain the condition at which 



pN +AN 

AE = Sfi(N)dN = 0. 

Jn 



(3.6) 



This condition is satisfied e.g. when the external chem- 
ical potential reaches the value corresponding to point 
1 in the inset of Figure 0b. At this point the integral 
of the difference between the dashed line, representing 
\x (A), and the bold line, corresponding to aeV g is equal 
to zero. This condition is similar to the Maxwell's rule 
well-known in the the theory of real gasesEl. Hence, in- 
stead of entering one by one (following the curve 1-2-3-4) 
a few electrons tunnel to the island simultaneously, in 
bunch. This bunching is brought about by attraction be- 
tween the electrons that start forming a new shell. As the 
hard disk model is particle-hole symmetric the bunching 
is also reproduced when electrons are about to finish the 
shell, as it is shown in the lower part of the inset. 

Comparison of Fig. |r|b with Fig. [| and [| shows that our 
model somewhat overestimates the oscillations of E and 
\x. For example our model predicts that SE oc N, while in 
the numerical experiments SE oc A 7 , with 7 = 0.8 ± 0.1. 
This overcstimation is associated with the main assump- 
tion of this model about the existence of the perfect disk- 
like shells, that brings about a strong degeneracy of the 
addition one-electron energies. In other words the one- 
electron density of states in our model comprises the se- 
ries of sharp 5-function like peaks. This is a good ap- 
proximation for small N, when these shells indeed exist, 
however in the form of perfect hexagons. At large N 
(N > 55) the configurations with perfect shells are never 
formed and the surface of the island is always rough. The 
peaks in the density of states existing at small N become 
smeared at large N. Nevertheless the density of states 
retains periodic variations. Hence while the variations 
of the ground state energy persist their amplitude is re- 
duced, with the exponent 7 going down from 1 to w 0.8. 
The detailed study of these variations is a subject of fu- 
ture work. 
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IV. THE COULOMB INTERACTION 



V. CONCLUSIONS 



The universality of the observed fluctuations of energy 
is not accidental. The case of the long-range interac- 
tions can be reduced to the simple model described in 
the previous Section. To do so one has to notice that 
if the number of electrons in the dot is not too big the 
electrons are situated in an effective potential that is very 
close to parabolic. This effective potential is significantly 
reduced with respect to the bare one due the to screening 
by the Wigner crystal: 



A^A~A/\e(l/R)\, 



(4.1) 



where A is the renormalized external potential, e (g) is the 
static dielectric function rendered by the Wigner crystal. 
In the Thomas-Fermi approximation this dielectric func- 
tion is given by 



e (q) = 1 + U q {d£ COII /dn)~ 



(4.2) 



Here U q = 2ire 2 / nq is the Fourier transform of the bare 
interaction potential, £ corr {n) = —(3e 2 n 1 / 2 /n is-the cor- 
relation energy of the Wigner lattice per particle a, w ith (3 
being a numerical constant. Putting formulas (4.1) and 
(4.2) together we obtain the following estimate for the 
renormalized confinement: 



A 



R 



A 

N 



(4.3) 



where ao = y 2/no\/3 cx N 1 / 6 is the lattice spacing in 

the middle of the island, and R cx N 1 / 3 is its radius (see 
also Appendix |Aj). Substituting A into Eq. (3.4) and 



replacing d by the characteristic lattice spacing in the is- 
land ao it is possible to obtain the uniform upper bound 
for the fluctuations of the total energy in the Coulomb 



\5E(N)\ < 0.3(e 2 / K ) 2/3 A 1/3 N 1/6 . 



(4.4) 



The numerical results for SEJN) presented in Fig. [I] agree 
with the prediction of Eq. (4.4) in the sense that 5E(N) 
has a very weak dependence on N. However the ampli- 
tude of the oscillations observed in the numerical exper- 
iments is by a factor of 5 smaller than predicted by this 
upper bound (see Fig. |l|). This discrepancy is due to 
the presence of the elastic deformations in the Coulomb 
Wigner crystal. In contrast to the case of the short-range 
interactions due to the smallness of the Young's modulus 
of the Coulomb crystal the newly added electrons on the 
surface of the island (shown by the dark circles in Fig. [?]) 
sink into the island, further decreasing the shell effect. 
This effect warrants a more detailed study. 



We have studied the system of classical interacting par- 
ticles in the parabolic confinement. We have observed 
the periodic variations of the ground state energy that 
are universal, i.e. independent of the exact form of in- 
teractions. These variations are shown to be of purely 
geometric origin and are associated with the combina- 
tion of the shell effect and the polaronic effect based on 
the confining parabola. If the electron-electron interac- 
tions in the island are short-range these variations bring 
about periodic bunching in the charging spectrum. 
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APPENDIX A: SMOOTH PART OF THE 
ENERGY OF THE WIGNER CRYSTAL ISLAND 
WITH COULOMB INTERACTIONS. 

In this Appendix we obtain the scaling of different con- 
tributions to the total energy of a circular dot filled with 
the Wigner crystal. In essence we derive Eq. (2.5). An 
analogous problem was considered in literature, for the 
system of point charges of the surface of spherea. 

We start by naming the length scales involved. They 
are R cx iV 1 / 3 , the radius of the island given by Eq. (|2.4|), 
and the lattice spacing in the center of the island ao = 

^2/noV3 cx N-V 6 (see also Eq. (gg)). 

The expansion of the total energy into the series of 
powers of N can be written in the following form: 



E„ 



Ees + E c . 



Eos + E s - 



(Al) 



The various terms in this expression are respectively: the 
electrostatic, the correlation, the overscreening, and the 
surface energies. The electrostatic energy can be calcu- 
lated from the solution of the electrostatic problem (|2.4|): 



ES 



2-Krdrn (r) (0o + Ar 2 ) 



(A2) 



where <po = 2AR is the total potential inside the island. 
As a result we obtain 



E ES = \NAR 2 
5 



(A3) 



The correlation energy describes lowering of the en- 
ergy due to the transition from liquid, assumed in the 
electrostatic model, to crystal. In case of a big island the 
energy reduction can be calculated in the local density 
approximation from tha-,known correlation energy of the 
uniform Wigner crystalH 



G 



■"corr 



27rrdrn (r) /?e 2 \/n(r)/ n 



(A4) 



12/1/V / ARe 2 
7tt 



where the coefficient /? was calculated in Ref. |j| For 
the triangular lattice it is equal to /3a = 1-960517, for 
the square lattice /?□ = 1.949299. The actually observed 
value of this coefficient j3 — 1.956 lies between these two 
values. 

n(r) 
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FIG. 8. The density profile in a small crystalline sam- 
ple (solid line). It is quite different from the solution of 
the electrostatic problem (dashed line) valid in the limit 
of very large island. The real density looks like a box if 
the sample contains only few crystalline rows in it. As an 
example on the bottom we show the ground state config- 
uration of N = 45 electrons. 

The next term in the energy expansion takes into ac- 
count the deviation of the density profile of the crys- 
talline island from the semicircle (2.4) obtained in the 
electrostatic model. Due to the screening the variations 
of the average density of crystal are slightly bigger than 
obtained in electrostatics. Note that in case of electron 
liquid these variations are weaker than the electrostatic 
ones. To emphasize this difference between liquid and 
crystal we call the corresponding correction to energy 
the overscreening energy. To estimate the correction we 
notice that the average density of the crystal deviates 
from the electrostatic density by Sn/n ~ —r s /R, where 
r s is the screening radius of the Wigner crystal equal to 
the negative lattice spacing: r s ~ — ao (Ref. As a re- 
sult the distances between electrons change by the same 
factor and the correlation energy acquires the following 
correction: 

Eos ~ E COII a /R ~ (e 2 / 'nf 3 A^N 2 '\ (A5) 

The last contribution to the total energy is the surface 
energy. It describes the deviations of the correlation en- 
ergy of electrons on the surface from that in the bulk. 
To find this deviation it is essential to know what is the 



lattice spacing at the edge of the island A (see Fig. ||). 
Below we address this problem in some detail. 

The electrostatic formula (^4|) is correct on the scales 
significantly exceeding the screening radius. As it was 
pointed out before the screening radius of-.the Wigner 
crystal is of the order of its lattice constantEEl. Hence the 
electrostatic formula breaks down at the distances from 
the edge 8R of the order of the lattice spacing on the 
surface (see Fig. 0): 



8R ~ A. 



(A6) 



Matching the electrostatic density of electrons at this dis- 
tance with the density of the Wigner crystal at the edge 
we obtain: 



n {R — A) = no 




R A 2 ' 



and 



A~^ /E \R 1/5 «^" 1/15 - 



(A7) 



(A8) 



Our numerical data agree very well with this theory and 
give the following coefficients: 



A = 0.8800' J R L 



5R = 0.764 /5 R}/ 5 . 



(A9) 



Notice that the inter-electron density at the edge differs 
very slightly from the lattice constant in the center of 
the sample if radius of the island is not too big. Hence 
the concentration of electrons in this case also does not 
change too much. This is displayed by a box like density 
profile in Fig. ||. 

The surface energy can now easily be estimated. It is 
of the order of the correlation energy per particle on the 
surface times the total number of particles there R/X: 



^_R 

k\ X 



(e 2 / K ) 2/3 A^N 7 ' 1 



(A10) 



Putting Equations (|I|) . . . flAlOp together we obtain 
the expansion (|2.5|). 
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